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Abstract 

The problem of jtt — ► 7T7T is studied using the axial anomaly, elastic unitarity, 
analyticity and crossing symmetry. Single variable dispersion relation is assumed. 
Using elastic unitarity relation, an integral equation equation for the lowest partial 
wave amplitude is obtained. The solution for this integral equation is obtained 
by an iteration procedure corresponding to that obtained from the vector meson 
dominance model but with multiple scattering effects taken into account. 



*Talk given at the "Fourth Workshop on Quantum Chromodynamics" 1-6 June 1998, the American 
University of Paris, Paris, France, 
t unite propre 014 du CNRS 



1 



1 Introduction 



One of the fundamental calculation in particle theory is the 7r° — > 77 decay rate It is a 
combination of Partial Conserved Axial Current (PCAC) and the short distance behavior 
of Quantum Chromodynamics (QCD): 

with 

e 2 iV 

F 77 = TT^f = 0.025GW- 1 (2) 

where e is the electric charge, f n = 93MeV and iV c = 3 is the number of color in QCD. 
This calculation is valid in a world where the 7r° is massless. Some corrections have to 
be made in order to take into account of the finite value of the pion mass. It turns out 
that the massless pion anomaly formula is in very good agreement with the pion life time 
data, impying that the correction to the physical pion mass is very small. 

Another Axial Anomaly result is the process 'jtt — ► tttt or its analytical continuation 
7 — > 37r 0. This last process requires more corrections because, for practical consider- 
ation, measurements will be done at energy far from the chiral limit where the anomaly 
formula is applicable; furthermore the analytical continuation from one process to the 
other is a delicate procedure due to the presence of the complex singularity which is ab- 
sent in the former reaction. The calculation of the process 77? — > ttti is in itself interesting 
because of future experiments being proposed at various accelerator facilities and also of 
its important role in the calculation of 7r° — > 77* /citeTruongl. 

The 77r — > 7T7T amplitude is given as: 

A( 7 (A;)7r> ) -> 7r + (pi)7r-(p 2 )) = iF 3w (s, t, u)e^ T e^p 0u p la p 2r (3) 

The kinematics are as follows: s = (k + Po) 2 ,t = (k — pi) 2 ,u = (k — p 2 ) 2 with k refers 
to the photon 4-momentum and po refers to the neutral pion 4-momentum and p\ and 
P2 are those of the charge pions. In terms of s and the c.o.m. scattering angle 8, we 
have t = a(s) — b(s)cos8, u = a(s) + b(s)cosd with a(s) = (3m 2 — s)/2 and b(s) = 

l/2(s-ml)/(l-Aml/s). 

In the chiral limit (the zero limit of the pion 4- momenta), the matrix element is given 
by the anomly equation: 

F ^(°) = X = t4t¥ = 9.7GW- 3 (4) 

where the zero in the argument of F 3n refers to the chiral limit of the massless pions; the 
number of colors iV c is equal to 3. 

Experimentally, A is measured at an average photon pion energy of OAGeV and as- 
suming that there is no momentum dependence F 37T , it is equal to [f|: 

X expt = 12.9 ±0.9 ±0.5GeV- 3 (5) 
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The agreement between experiment and theory is not very good. Taking into account 
of the momentum dependence of F 37r , the following value of A is obtained ||: 

X expt = n g ± q g ± Q 5Ge y-3 (6) 

There is still disagreement between theory and experiment. 

The calculations of this process are usually done within the Vector Meson Dominance 
models (VMD) [|7|, || |§; recently it is discussed within the framework of Chiral Perturba- 
tion Theory (ChPT) [TIJ and also a combination of ChPT and VMD ||. 



2 Integral Equation Approach using Elastic Unitar- 
ity Relation 

In this talk, the process 'jix — > ixix is studied using dispersion relation and elastic unitar- 
ity An integral equation, similar to the pion form factor Muskhelishvilli-Omnes integral 
equation |12j is obtained. The difference is that the integral equation to be treated here 
is much more complicated due to crossing symmetry; no exact solution has been found. 
We shall get the solution of this integral equation by an iterative procedure, but with the 
crucial property that the iterative solution satisfies the phase theorem at every steps as 
required by unitarity. 

We make here the same assumption as in the pion form factor calculation. Because 
the pion form factor can be calculated to within 15 % (in the amplitude) Jill, we ex P ec t 
to have the same degree of reliability for this new integral equation. 

If a better accuracy is demanded, similar to the calculation of the pion form form 
factor, a more precise low energy measurement for '-fix — > ixix is needed. 

The strong P-wave ixix scattering phase shifts are supposed to be known and are 
given by the experimental data up to 1 GeV or higher which show the existence of the 
p resonance at 0.77 GeV with a width of 0.151 GeV. We shall not directly make an 
assumption on the Vector Meson Dominance, but try to find a solution consistent with 
the constraints of the elastic unitarity, crossing symmetry and also of the low energy 
theorem, Eq. (f|). 

The process '-fix — > ixix is a completely symmetric reaction in the three variables s,t,u 
i.e. the same amplitude describes not only the reaction '-fix — > ix + ix~, but also the two 
other amplitudes involving the permutations of the pions. It is assumed that the scatter- 
ing amplitude -fix — > ix + ix~ can be represented by a single spectral function dispersion 
relation: 

^/ x t- A s — m l) f°° a(z)dz , , , , , 

F(s,t,u =A+" ~ ~, £f -} + [s~t] + [s~u] (7) 

IX J Ami \ z ~ m i)\ z — s ~ le ) 



where for simplicity we drop the subscript 3ix in F(s,t,u) and set A = F(m^,m^,m^), 
the value of the scattering amplitude at the symmetry point; the relation between A and 
A will be discussed below. Eq. (0) shows explicitly the symmetry in s, t, u variables. 
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Because the p resonance occurs at a much higher energy than m%, we expect that 
A differs little from A. We can first estimate the value of A by examining the large N c 
limit of QCD. In this limit, the pion loop corrections are neglected and the vector meson 
propagator appears as a pole with zero width. Introducing a contact term to satisfy the 
anomaly equation, the Vector Meson model yields A = A(l + 3m^/2s p ) ~ 1.05A. We shall 
use this value in the following. 



2.1 Partial Wave Amplitude and Integral Equation 

Taking into account of the lowest partial wave projection of F(s, t, u), F(s), together with 
the strong P-wave nn — > nn amplitude, the elastic unitarity relation gives: 

a(s) = F{s)e- i5{s) stn5{s) (8) 

where 5 is the P-wave irn phase shift obtained from the available experimental data which 
show that they pass through 90° at the p mass with a width of 151 MeV and that there 
is no measurable inelastic effect below 1 GeV. We have 



F( 8 ) = a, ( s " m D r F{z)e~ iS{z) sin 6(z) ^ 



7T J Ami ( z ~ m 1)( Z — S — ie) 

1 f 00 v -is( z ) . cv u 1 , i z - a(s) + b(s) . 2 
+- / F(z)e l5{z > sm5(z){—ln ^ — — -}dz 

TT J Ami K S ) Z — a(s)—b(s) [Z — 171%) 

(9) 

Eq. (§) is a complicated integral equation. It is similar to, but more complicated 
than the Muskelishvili-Omnes (MO) type [I2[, because the t and u channel contributions 



are also expressed in terms of the unknown function F(s). It should be noticed that the 
first term has a cut from 4mjj. to oo and the second one has a cut from to — oo. This 
remark enables one to solve the integral equation by the following iteration scheme which 
converges very fast. 

The iterative and final solutions can be expressed in terms of the function D(s,0), 
normalized to unity at s = and defined in terms of the phase shift 5 as given by: 

i s p 6{z)dz 

exp— / — — (10) 



D(s) 7r JAmi z{z — s — ie) 

Other functions D normalised to unity at s = Sq can be expressed in terms of the function 
D(s,0) by the simple relation D(s, s ) = D(s,0)/D(sq,0). 



2.2 Iterative Solutions 

As it was remarked above, the Integral Equation has both right and left cuts. Because 
of this analytic structure, we can define an iteration procedure which consists in splitting 
Eq. (|9|) into two separate equations: 
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f« w = i + Tt%) + ^ r ry* 1 "! * (id 

3 7r J4m| (2; — m~)[z — s — ze) 



and 

Tt%) = f + " T F^ We ^.m* W { * to | f = a ' 3 > + j£> | -—^U, 
3 7r J4m| b(s) z — a(s) — b(s) {z — m%) 

(12) 

where % > 1 and is the value of the function F(s) calculated at the i th step in the 
iteration procedure; the Born term T B ~ (s) is calculated at the i th — 1 step. The Born 
term is real for s > and has a left cut in s for s < 0. In writing Eqs. (|Tli|l2l), care was 
taken to preserve the symmetry in the s, t, u variables for the function F(s, t, u) which 
requires us to split symmetrically the subtraction constant A in Eq. @ into three equal 
pieces, one contributes to Eq. (|IID the other two to Eq. (|Tj). 

The solution of the integral equation Eq. (|TT]) is of the MO type [fT2 |: 



U 3D( S ,m2) + B U + £>0,O 7T Am' ( z _ m 2)( 2 _ s _, e ) 

(13) 

where it is assumed that the well-known polynomial ambiguity inherited in the MO inte- 



gral equation is absent except the constant term which is required by the low energy 
theorem. It can be shown that the phase theorem for (s) is satisfied. The first term on 
the R.H.S. of Eq. (|13|) represents the p vector meson contribution to F^'(s), the second 
term is roughly the same as the ptfy vertex correction due to the exchange of a resonant 
pair of P-wave pions (p vector meson). (The higher polynomial ambiguity could represent 
some uncalculable inelastic effect occuring above the inelastic threshold). 

One arbitrarily defines the convergence of the iteration scheme at the ith iteration 
step when | | / | F^ 1 ' | differs from 1 by less than 1% or so in the energy range 
from the two pion threshold to lGeV. (Alternatively one can also require that the ratio 
| T B | / | Tg I to be unity within an accuracy of 1% or so). 

Once the solution for the partial wave is obtained we should return to the calculation 
of the full amplitude. This can be done by combining the Tg 1 ' Born term in Eq. (|T3| ) 
with higher uncorrected partial waves (for rescattering) from the t and u channels to get 
the final solution: 

F»( a , t, u) = ^[{^^(1 + S/^ 1 ^))} + {{s - t)} + {(s - «)}] (14) 
where the function denotes the multiple rescattering correction: 

7r J 4mJ (z — ml)(z — s — ie) 



Projecting out the I = partial wave from Eq. fllil) , we arrive at Eq. (]T3|) with T l B (s) 
replaced by T B (s). Because of the assumed criteria for the convergence of the iteration 
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scheme, T]j _1 (s) ~ T % B (s) it is easily seen that FW(s) has the phase 5, using the result of 
Eq. (0). The remaining (even) higher partial waves I > are all real because we have 
assumed that the strong final state interaction of the higher partial waves are negligible. 
The final solution Eq. ([TJD is completely symmetric in the s, t, u variables. 



3 Numerical Solutions 

In order to carry out the iteration procedure to find the solution of the integral equation 
one has to parametrize the function D(s,0), normalized to be unity at s — 0, in terms 
of the experimental P-wave phase shift. In the following we shall make three different 
parametrisations for the D function. They are all expressed in terms of the P-wave hit 
phase shifts 5 as given by Eq. (|10|). 

The function D _1 (s) can be parametrised as follows fT4|, ITT]: 



D (S) l-s/s Rl -^{(s-4ml)H wn (s) + 2s/3} (16) 

The p mass is defined as the vanishing of the real part of the denominator and is equal 
to sr in the narrow width approximation. Using ^/s^ = 0.770(jeV in Eq. fllBD, we 
have T p = 155.6MeV for /„. = 0.093GeV . which is very near to the experimental value 
of 150.7 ± 1.2MeV. If we want to fit the -Di(s) to the experimental width, we can 
phenomenologically change f n = 0. 0945GeV in Eq. fll6|) . 

In the above subsection we have discussed the comparison between our result and the 
low energy experimental data. Our result is about 1 standard deviation too low if the 
the value of A = 12.9 ± 0.9 ± 0.5GeV~ 3 is used. Taking into account of the momentum 
dependence of F 3lT then the experimental value of A is 11.9 ± 0.9 ± O.SGeV^ -3 and is 
in excellent agreement with our calculation. At low energy, the effect of the multiple 
scattering is important and the good agreement with the experimental data is obtained 
thanks to this effect which is fully taken into account in solving the integral equation. 

We want to discuss now the result of our integral equation in the p resonance region 
where there are some experimental measurements. The experimental data are usually 
analyzed in term of the p — > 7T7 width using the Breit-Wigner formulae: 

«(^^n\ S p T( P -> 27T)r(p -> 7T 7 ) 

<T(77T -> p) = 2 r^- -= (17) 

where T p is the total width of the vector meson p. In terms of the matrix element for 
77r — ► 7T7T, F 37r (s,t,u), we have: 



(77r — > titi) = — (s — m w ) sin 9 (1 



dcos9 y ' 1024tt 

In arriving at Eq. ( fI7|) we must assume that there is no background term in the cross 
section which might interfere with the Breit-Wigner term which may not negligible. By 
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background term we mean amplitudes which do not have the Breit Wigner form which 
might be present. In our formulation the background term from the t and u channels 
are automatically taken into account as can be seen from Eq. flT8"|). It turns out because 
the p width is fairly small, the Breit- Wigner approximation, Eq. (|17D , is accurate to few 
percents, and hence we can use it with confidence. Our method is therefore to equate Eq. 



(Jig) and Eq. ( [171) to calculate the p — > 7T7 width which is a convenient way to express 
our solution in terms of physically measurable quantity. We shall use this method when 
we compare our result with those given by other models. 

Using the D-function as given by Eq. ( |T0"D which satisfies the KSRF relation, we 
obtain: 

T(p -> 7r 7 ) = Q3KeV (19) 

This value is to be compared with the world average for the charge p — > 7r 7 width , 
T(p + — > 7r + 7) = 68 ± 7KeV. The corresponding for the the neutral p is very different, 
r(p° — > 7r°7) = 120 ± 30KeV. As explained above, there should be no difference between 
the values of the charge and neutral radiative widths. 

Using T(lj — > 7r 7 ) = 716 ± 42KeV and the SU(3) relation, we get T(p — > ttj) = 
80 ± 5KeV which is about 20% larger than the present data. 

How accurate is our calculation with the assumption of the elastic unitarity relation? 
From our experience of using the same assumption for the pion form factor calculation, 
the maximum value of the square of the absolute value of the pion form factor at the p 
peak is 32 which is 20-25% too low compared with the data. We expect to commit the 
same magnitude of error using the elastic unitarity for our problem. If the experimental 
data on p — > itj width was changed to the SU(3) value of 80KeV, our value would still 
be consistent with the data. 

To put it more quantitatively, the pion form factor V(s) can be mulptiplied by a real 
polynomial, say (1 + 0.13s/s p ), which phenomenologically represents the inelastic effect 
occuring at a higher energy. This modification gives an excellent fit to the experimental 
data on the time- like pion form factor up to lGeV 2 . 

We can likewise multiply our solution by the same factor which might represent some 
correction to our assumption of the elastic unitarity relation for our problem. This dis- 
cussion is purely phenomenology and is not needed at present because our calculation is 
in good agreement with data. More accurate measurement of F 3n at low energy will be 
necessary in order to limit the size of this term. 

In conclusion,the process 7 7r — > tttt is calculated using the low energy theorem, ana- 
lyticity, elastic unitarity and crossing symmetry which are fundamental conditions for a 
theory involving strong interaction. The final result shows that one effectively takes into 
account of the (unstable) p model in the s, t and u channels and their rescattering effect 
treated in a self-consistent way. It is important to put these results to experimental tests. 
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